Floquet Spectrum and Transport Through an Irradiated Graphene Ribbon 
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Graphene subject to a spatially uniform, circularly-polarized electric field supports a Floquet 
spectrum with properties akin to those of a topological insulator, including non-vanishing Chern 
numbers associated with bulk bands and current-carrying edge states. Transport properties of this 
system however are complicated by the non-equilibrium occupations of the Floquet states. We 
address this by considering transport in a two-terminal ribbon geometry for which the leads have 
well-defined chemical potentials, with an irradiated central scattering region. We demonstrate the 
presence of edge states, which for infinite mass boundary conditions may be associated with only 
one of the two valleys. At low frequencies, the bulk DC conductivity near zero energy is shown to be 
dominated by a series of states with very narrow anticrossings, leading to super-diffusive behavior. 
For very long ribbons, a ballistic regime emerges in which edge state transport dominates. 

PACS numbers: 72.80.Vp,73.23.-b,73.22.Pr 
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Introduction and Key Results - The electronic proper- 
ties of graphene are very unusual among two-dimensional 
conducting systems, in large part because the low energy 
physics is controlled by two Dirac points, which form the 
Fermi surface of the system when undoped [iHl] . One of 
the very interesting possibilities for this system is that, 
with spin-orbit coupling, it may represent the simplest 
example of a topological insulator Topological 
insulators are systems for which the bulk spectrum is 
gapped, but which support robust, gapless edge states. 
Unfortunately, spin-orbit coupling in graphene appears 
to be too weak to allow observation of this behavior with 
currently available samples. 

Very recently, theoretical studies have suggested that 
an analog of topological insulating behavior can be in- 
duced in g raphene by a time-dependent electric poten- 
tial |7H12|. The proposal entails exposing graphene to 
circularly-polarized electromagnetic radiation of wave- 
length much larger than the physical sample size, such 
that only the electric field has significant coupling to the 
electron degrees of freedom. The periodic nature of the 
field necessitates that the quantum states of the elec- 
trons are solutions of a Floquet problem, characterized 
by a "quasi-energy" Eq. with allowed values in the inter- 
val [— (x'/2, aj/2], where a; is the frequency of the radiation 
[isj . This same physics allows one to induce topological- 
insulating properties in a variety of systems that are oth- 
erwise only "almost" topological insulators jsj. Because 
time-reversal symmetry is explicitly broken in this sys- 
tem, it should support a Hall effect Q which, when mea- 
sured in an appropriate geometry, may be quantized . 

A key challenge one faces in determining transport 
properties of this system is the assignment of electron 
occupations to the Floquet states. In general, the 
quasienergies Eq, cannot be simply inserted as energies 
in a Fermi-Dirac distribution since they are limited to a 
finite interval of real values, determined by the frequency. 
In this study, we assume Fermi-Dirac distributions only 
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1: {Color online) Schematic diagram of device geome- 



for the incoming waves far in the non-irradiated leads. 
Thus we assume that electrons are injected and removed 
from the system via highly doped, ideal leads in which 
any possible effects of an electric field have been screened 
out, and that the transport within the (finite) irradiated 
region is quantum coherent. Our geometry is a direct 
analog of one studied in Ref. [13 for the time-independent 
case, and is illustrated in Fig. [TJ Identical leads are taken 
to be made of highly doped graphene. This leads to a 
vanishing time-averaged current in the absence of a DC 
bias, which avoids photovoltaic (charge pumping) effects. 

For an infinite ribbon geometry, the momentum kx 
along the ribbon axis is a good quantum number, and 
one may compute an effective band structure for the Flo- 
quet eigenvalue Fig. [5] illustrates this for several 
cases. Fig. [2][a) displays the bands closest to e = 
for a relatively large frequency, hw — it where t is the 
tight-binding hopping parameter. This result is illustra- 
tive for its relative simplicity, but is only relevant to very 
small system sizes for which the electric field may be ap- 
proximated as uniform throughout the sample [111] . One 
may see two sets of minima/maxima, corresponding to 
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FIG. 2: [Color online) (a) Tight-binding Floquet spectrum 
as a function of for a graphene ribbon oriented in the 
zigzag direction, subject to a circularly polarized electric field 
of magnitude Eq with eEoa/hw = 0.5, hw/t = 3.0, where a 
is the Bravais lattice constant and t the hopping parameter. 
(b,c) Spectra from Dirac equation with infinite mass bound- 
ary conditions, corresponding to two different valleys. Only 
one valley carries edge states. Frequency iua/vp = 5, electric 
field amplitude Eo/hw = 1, ribbon width W = 45a. 



the two valleys, separated by a gap that does not vanish 
even as the system width W becomes very large. This 
intrinsic gap in the Floquet spectrum is an analog of that 
which opens in the presence of spin-orbit coupling for the 
static graphene system In further analogy to this, a 
pair of edge states traverses the gap and connects the two 
valleys, while (anti-)crossing around — 0. 

Figs. IH^bjC) illustrate corresponding results for Flo- 
quet spectra of two individual Dirac points subject to the 
circulating electric field, corresponding to the two differ- 
ent valleys, with infinite mass boundary conditions. In 
analogy with a topological insulator, one sees that both 
valleys develop gaps which do not go away as W becomes 
large, but only one supports edge states, so that the num- 
ber of states at an edge is unaffected (modulo 2) by the 
change of boundary condition, while the details of how 
these states disperse are changed. Since the system with 
its edge state may be well-described qualitatively with 
just the single valley illustrated in Fig. [2][b), we focus 
our attention on this case for the transport calculations. 

Fig. [3] illustrates typical results for the conductance of 
the system as a function of the scattering region length, 
L. At large L one can see the conductance level off to 
a constant value, indicating the presence of edge states. 
For sufficiently wide samples this ballistic behavior will 
be robust against disorder, since states carrying current 
in opposite directions reside on opposite sides of the sam- 
ple. The presence of a finite current for large L is in 
marked contrast to the behavior in the absence of the 



+ "'I, 



G{e^/h) 




0.1 



L/a 



FIG. 3: (Color online) Conductance vs. L for several differ- 
ent bandwidths: Qn^^^W = vrrimax, with rimax ~ 10 (blue as- 
terisks), 15 (lavender open squares), 25 (aqua closed squares). 
For these plots, W — 20a, uja/vp = 5, and number of time 
steps is 13. Note approximate power law behavior G ~ 
with b « 0.65 for L < W. For comparison, results for Eq — 
displayed as red crosses, displaying G ~ 1/L behavior for 
L < W . Inset: Illustration of Sa vs. for infinite system 
with periodic boundary conditions, showing Floquet copies of 
spectrum and resulting avoided crossings for kx 7^ 0. 



radiation, for which the conductance vanishes [T4|. Sur- 
prisingly, the conductance exceeds the unirradiated con- 
ductance for all values of L, in spite of the fact that a 
gap has opened in the (Floquet) spectrum of the bulk 
system Q • The explanation of this lies in the fact that 
Floquet eigenvalues are restricted to a finite interval, for 
example — w/2 < < uj/2, so that energy states out- 
side this interval in the absence of the periodic potential 
are folded into it. For fixed ky, this results in a series 
of repeated crossings at values of v-pk^ ~ muj, where m 
are non-zero integers. In a system without edges (i.e., if 
one considers periodic boundary conditions rather than a 
ribbon) these become avoided crossings with very small 
gaps, as we explain below, which transport current via 
evanescent states. Remarkably, transport through these 
states results in super-diffusive behavior, G ~ l/L'^, with 
b < 1, as is apparent in Fig. [31 This non-analytic be- 
havior reflects an explosive growth of the decay length of 
evanescent states with large |to|, which has the form 



|r7i| U! 



^l^^'l —ri\m\ 



\m\ uj ' 



(1) 



where Aq = Eq/lj characterizes the electric field ampli- 
tude, assumed to be small, and 77 is a number of order 
unity. As we argue below, the rapid growth of with m 
results in an anomalously large penetration of the elec- 
trons into the irradiated region. 

Transmission Through Irradiated Region - The wave- 
functions for Dirac electrons subject to circularly po- 
larized radiation obey the time-dependent Schrodinger 
equation [—idt + H]'i' = Hp'i' = 0, where the Hamilto- 
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nian of the system has the form 
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Px 'f'Py H~ -^x '^^y 







(2) 
(Note 
1 in this ex- 



where A = Ao(cosuj<, sinwi), and Px.y 
we have set the Fermi velocity up 
pression.) Since H is periodic in time, the solutions 
will be Floquet states, which have the form ^(r,t) — 
e^^°*[$A(r,t),$i3(r,t)]T, with $^(r,i+T) = $^(r,t) and 
T = 2t:/uj. Adopting infinite mass boundary conditions 
leads to the conditions '14] $^(a;, y — 0,t) = ^B{x,y — 
0,t) and ^Aix,y = W,t) = -$s(a;, y = W,t). Eigen- 
states of H which meet these boundary conditions are 



<T,(n,s) — 



sz 

1 - 



' SZ J \ sz — 1 

^ik^x — iAyy 



(3) 



where z = {qx + iqn)/q, = - A^, g„ = (n+i)7r/W^, 
s = ±1, H$ = sg$, and Af^ is a normalization constant. 
States in the leads of the system are generated by setting 
A = in Eq. El 

Our strategy for finding the Floquet eigenvalues is 
to discretize time and expand the Floquet operator 
Hp in instantaneous eigenstates of the Hamiltonian T-L. 
Writing = (i) Jt.^t, we may then write 

{ni,Si,ti\'H{t)\n2,S2,t2) = E(^ni,st)itl) Sni,n2Ssi,s-2^ti,t2- 

The Floquet operator can them be written as a matrix 
of the form 



i 

~2At 



(ni,si,ti|n2,S2,ti - At) 6t2„ti-At 
-{ni,si,ti\n2,S2,ti + At) St2,ti+At , (4) 



whose eigenvalues are the allowed values of Sa for an infi- 
nite ribbon. Note that the states are implicitly functions 
of kx ■ Diagonalization of Eq. |4] generates results such as 
those depicted in Fig. [DJb). These results were obtained 
for < n < 30 and 19 time slices. 

Turning to the conductance, in the leads there are no 
microwaves, so t hat eigen values of Hp with A = have 
the form e = ±\/kl + ql + smimuj At) / At = ±En{kx) + 
e^. This implicitly defines an equation for k^^ which 
depends on the integers n and m. In the scattering region 
as well, for a given subband of the Floquet spectrum, we 
need to know the values of k^ that will give some specified 
Floquet eigenvalue £„. This is equivalent to finding the 
values of k^ where the bands illustrated in Fig. [^fb) 
cross some specified horizontal line. Note that if there is 
no such crossing for a given band then the corresponding 
kx is actually complex, indicating an evanescent state. 

In order to match the wavefunctions in the leads to the 
scattering region we need to know these latter values of 
kx , for a given Floquet eigenvalue . To accomplish this 
we multiply the eigenvalue equation by ax to obtain 



This is a non-Hermitian matrix equation which we ap- 
proximately solve in a manner analogous to what we did 
for the original Floquet equation. Eigenvectors give us 
the wavefunctions in the scattering region, and the eigen- 
values kx that enter into the plane wave part of the wave- 
function e**^^^. A comparison of these solutions to £a{kx) 
obtained by diagonalizing iJp reveals excellent agreement 
between the two calculations. 

We now have analytic formulas for the wavefunc- 
tions in the leads, and approximate numerical solutions 
for them, represented in a finite basis of the states 
<j)(n,s)(y) ^ik^(,n,m,s)x+tmt ^ -j^ Scattering region. These 

need to be matched at two junctions. To accomplish this 
we match the wavefunctions on discrete points in y, tak- 
ing Vj = (j + 5)VF/('^max + l), where rimax + l is the num- 
ber of transverse states retained, and j £ {0, . . . , n„iax}- 
This defines a set of linear equations which we solve nu- 
merically, and from which the matrix T^p^{E,E + e^), 
representing the time-averaged transmission across the 
structure, may be obtained. (Here q,p represent trans- 
verse channels in the left and right leads, respectively, 
and E is the energy of an impinging electron from the 
left.) One may show that by matching on these par- 
ticular points one enforces current conservation in the 
solution; this is confirmed numerically to 1 part in 10^. 
The conductance is finally given by [isl , [l6j 



(6) 



p,q,n 



[{-idt - ea)cTx + ia-zPy]i> = kxtp- 



(5) 



where Ep is the Fermi energy in the leads. 

Evanescent Transmission in Irradiated System - A 
prominent result from the calculations of the two- 
terminal conductance is an approximate power law be- 
havior G ~ when L < W. As b is non-integral 
this represents non-analytic behavior, and the fact that it 
emerges well above the large L value of G suggests it is a 
result of evanescent state transport. This behavior turns 
out to be rather natural when one accounts for higher 
order crossings at non-zero kx of the Floquet spectrum. 
In a realistic situation, for example if the impinging radi- 
ation is in the microwave regime, hu; ~ 10~^ eV, one will 
have many such crossings since the graphene bandwidth 
1 eV) is relatively large. 

To demonstrate the behavior, we consider a simpler 
problem in which there is a half-space of irradiated 
graphene, and a half-space of highly-doped, unirradiated 
graphene, joined across x — 0, and we adopt periodic 
boundary conditions in the transverse direction so that 
there are no edge states. For this situation ky is a good 
quantum number. States with zero energy are perfectly 
backscattered in this case because there are no propa- 
gating states with zero Floquet eigenvalue. In the steady 
state situation there is a charge density tail penetrating 
the irradiated graphene, which has an approximate power 
law falloff with x. 

To see this, we need to know how evanescent wavefunc- 
tions fall off with X inside the irradiated graphene. In the 
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limit of vanishing microwave amplitude the Floquet 
spectrum will have pairs of states crossing Eq, = at 



kx = ±fcm = ±^ {niujY — ky. The degeneracy is lifted 

for non- vanishing Aq, which we take to be small com- 
pared to uj. The resulting anti-crossing will occur at high 
order in perturbation theory, since the degeneracy occurs 
for states with time dependence exp(±iTOwt), whereas 
the perturbation V — Aq (ct^ cos ojt + sin tot) connects 
states whose frequencies differ by a single unit of uj. Thus 
the gap that opens at the crossing will be proportional 
to (^o/w)^™. Since the inverse gap is essentially the lo- 
calization length we wish to compute this. One approach 
is to use the resolvent operator G{z) = (z — H-p)^^ with 
Hp ^ Ho + V. Poles of the 2x2 matrix 



G = 



{m, 



G I m , 



-) ("^,-1^1 
-) {-m,+ \G\ 



m , +; 
- m ,+) 



(7) 



where s = ± indicates a particle-like or hole-like 
(-) state, are the eigenvalues of Hp. Expanding G in 
powers of V allows one to define a self-energy, G^^ = 
Gq^ — S = z — Hq — S, with Ho = for the states 
of interest. The diagonal components of S simply shift 
the precise location of the anticrossing on the k^ axis 
and may be ignored. To lowest non-trivial order, the 
off-diagonal components are I]± = {m , — | E | — m, -f ) = 
(— m , -f S TO, — )*, with 



S+ = 



{kx + iky)Ao 



2k 



2m 



(8) 



m—1 

X Jl \{n 

n=— m+1 



TL—l 

n [(n,-l-|G("'|n,+)-(n,-|G(°)|n,-) 



Evaluating the matrix elements and setting z — 0, one 

-1 



finds 



^ - 1 



In the 



|E±| = A2™(TOC^)l-2mJ|^^ 

limit of large to, the product can be evaluated; noting 
that = 2|S±|, one arrives at the estimate in Eq. [1] 
with f] = 4(1 - In 2) PS 1.227. 

To see the connection with power law behavior, one 
needs to develop matching conditions at the a; = inter- 
face with these wavefunctions. This is a tedious but in 
principle straightforward exercise [l7j , yielding the result 



that the contribution to the density from large n has the 
time-averaged form 



p{x) ^ J2 



2n 



where r is of order unity, and depends on ky. The sum 
may be estimated by assuming it is dominated by a single 
term at large n when x is large; maximization yields 



ln(a;a;) 



G 



where G is independent of L, £„ is a correction of 
order ln[ln(ti;a;)] , and y = — 21n(^o/i^) ~ Using 
this term to estimate the sum yields the result p{x) ^ 

(^2;)21n(a/c^r)/(y+2£„) gj^^^g ^j^jg ^ ^g^^j^ fuuctiou of 

ky (through t), we see that summing over the transverse 
modes should result in an approximate power law density 
tail. Such a fall-off is expected to lead to similar behav- 
ior in the transmission. It is interesting to note that the 
actual value of the exponent is relatively insensitive to ui 
and Aq, since these enter only through logs; this is consis- 
tent with our results, for which the observed power tends 
to remain in the interval 0.6 < b < 0.75 over a variety of 
choices for uj and Aq. 

In summary, two-terminal transport through an un- 
doped graphene ribbon subject to a circularly rotating 
electric field has a conductance that reveals the unusual 
nature of the Floquet spectrum of this system. Evanes- 
cent transport through relatively short ribbons is super- 
diffusive due to a series of near crossings with very small 
gaps. At larger ribbon lengths, transport becomes bal- 
listic, revealing the presence of edge states which are a 
hallmark of the topological nature of the spectrum. 
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